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Suzaku Observations of Metal Distributions in the Intracluster Medium 

of the Centaurus Cluster 
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Abstract 

We report the first observations of metal distributions in the intracluster medium of the Centaurus 
cluster up to ~ 0.17 ri 80 with Suzaku. Radial profiles of the O, Mg, Si, S, Ar, Ca, and Fe were determined 
at the outer region of the cluster, and their variations appear to be similar to each other. Within the cool 
core region (r < 0.045 riso), all the metal distributions sharply increased toward the center. In the central 
region (r < 0.015 riso)) the abundances of Si, S, Ar, Ca, and Fe were 1.5-1.8 solar, while those of O and 
Mg were approximately 1 solar. The derived abundance ratios of O and Mg to Fe were slightly lower than 
those of a set of other clusters. In contrast, the calculated mass-to-light ratios (MLRs) for O, Mg, and Fe 
were larger than those of the other clusters. For the outer region of the cool core (r > 0.07 riso), all the 
abundances were almost constant at 0.5 solar. The derived MLRs were comparable to those of the other 
clusters. This suggests that the cD galaxy of the Centaurus cluster efficiently supplies more Fe than the 
other clusters. 
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1. Introduction 

The metal abundances of the intracluster medium 
(ICM) provide important information regarding the chem- 
ical history and evolution of clusters. The metals in the 
ICM were primarily produced by supernovae (SNc) in 
early-type galaxies (Arnaud et al. 1992; Renzini ct al. 
1993). Si and Fe are synthesized both in SNe la and 
SNe II, while O and Mg are synthesized predominantly 
in SNe II. O and Mg abundances, therefore, are crucial 
for understanding the history of massive stars in clusters 
of galaxies. 

Because the Fe-K lines are prominent in the spectra of 
the ICM, the Fe abundances of the ICM have been studied 
in detail. ASCA first showed Fe distribution in the ICM 
(Fukazawa et al. 1998; Fukazawa et al. 2000; Finoguenov 
et al. 2000; Finoguenov et al. 2001; Ezawa et al. 1997). 
De Grandi and Molendi (2002) derived Fe distribution 
and its gradient from Beppo-SAX observations. Recently 
XMM-Newton and Chandra observations have been used 
to study spatial distribution of Fe up to 0.3-0.4 r^o 
(Vikhlinin et al. 2005; Baldi et al. 2007; Maughan ct 
al. 2008; Leccardi & Molendi 2008; Matsushita 2011). 
Within the cool core regions, Fe abundance in clusters 
decreases sharply toward the outer region. Matsushita 
(2011) found that outside the cool cores, both cD and 
non-cD clusters have similar Fe abundance profiles, which 
are flatter than expected from numerical simulations. A 
simple explanation is that a significant fraction of Fe was 
synthesized in an early phase of cluster evolution. 

XMM-Newton provided the means to study the O and 
Mg abundances in the brightest cool cores of clusters and 



groups of galaxies (Finoguenov et al. 2002; Matsushita et 
al. 2003; Tamura et al. 2003; Matsushita et al. 2007b). 
However, a higher background level and a strong in- 
strumental Al line of the MOS detector onboard XMM- 
Newton can cause problems in measuring O and Mg 
abundances outside cool cores. The X-ray Imaging 
Spectrometer (XIS) instrument (Koyama et al. 2007) on- 
board Suzaku (Mitsuda et al. 2007) offered an improved 
line spread function below 1 keV coupled with a lower 
background level. With the Suzaku satellite, O and Mg 
abundances of the ICM outside the cool cores of sev- 
eral clusters and groups of galaxies were measured up 
to 0.2-0.3 riso (Matsushita et al. 2007a; Komiyama et 
al. 2009; Sato et al. 2007a; Sato ct al. 2008; Sato ct al. 
2009a; Sato ct al. 2009b; Sato ct al. 2010). Within the 
cool cores, the derived abundance patterns of O/Mg/Si/Fe 
are comparable to the solar ratio, obtained by adopting 
the solar abundance by Lodders (2003). Combining the 
Suzaku results with SNe nucleosynthesis models, Sato et 
al. (2007b) showed the number ratios of SNe II to la to 
be ~3.5, and that Fe was synthesized mainly by SNe la. 
A slight increase in O/Fe and Mg/Fc with increasing ra- 
dius was detected in Abcll 1060 (Sato et al. 2007a) and 
AWM 7 (Sato et al. 2008), although the error bars were 
reasonably large. MLR, the ratios of metal mass in the 
ICM to the total galaxy (stellar) luminosity, is important 
for studying the origin of metals, because metals are syn- 
thesized in stars. The integrated values of Fe mass-to- light 
ratio (IMLR) profiles of the kT =2-3 keV clusters, such 
as Abcll 262, and AWM 7 observed with Suzaku increased 
outward from 0.1 riso to 0.3-0.4 riso. These results im- 
ply that Fe in the ICM extends farther than stars in the 
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outer regions. These results also indicate early metal en- 
richment process in the ICM. 

The Centaurus cluster (Abcll 3526) is a nearby X- 
ray bright cluster with a cool core (Allen and Fabian 
1994; Ikebe et al. 1999; Sanders et al. 2008; Takahashi 
et al. 2009) at a redshift of z = 0.0104. At the clus- 
ter center, the cD galaxy NGC 4696 is associated with 
the low-power radio source PKS 1246-410 (Taylor et al. 
2006). The cluster has a plume- like structure around the 
cD galaxy. The cluster has a subcluster structure Cen 45 
at southeast of cluster center. With ASCA observations, 
the Cen 45 region has a hotter temperature than sur- 
rounding regions (Furusho et al. 2001). The Centaurus 
cluster has the highest central Fe abundance of the ICM 
of ^2 solar among nearby cool core clusters (De Grandi 
and Molcndi 2009; Matsushita 2011). From a Chandra ob- 
servation of the Centaurus cluster (Sanders et al. 2006), 
the abundances of O, Ne, Mg, Si, S, Ar, Ca, Fe, and Ni 
were studied up to 0.035 r 180 . With XMM- Newton, the 
abundances of O, Si, and Fe up to 0.1 r\so and abundance 
of Mg within 0.03 rigo were derived (Matsushita et al. 
2007b). The central abundances of O and Mg are a fac- 
tor of two lower than those of Si and Fe. These results 
indicate that Fe in the central region is mainly supplied 
by the ejecta of SN la in the cD galaxy. 

In this paper, we report the abundance pattern of 
O/Mg/Si/S/Ar/Ca/Fe of the Centaurus cluster observed 
with Suzaku up to 0.17 ri§o within and outside the cool 
core. We adopt the solar abundance values given by 
Loddcrs (2003), in which the solar Fe abundance relative 
to H is 2.95 x 10 -5 in number. We use the Hubble constant 
Hq = 70 km s _1 Mpc -1 . The luminosity distance to the 
Centaurus cluster is Z?l =44.9 Mpc, and V c orresponds to 
13 kpc. The virial radius r lso = 1.95/i r0 1 oV / (fcT)/10 kcV 
(Markevitch ct al. 1998; Evrard et al. 1996) is approxi- 
mately 1.74 Mpc for the average temperature of (kT) = 
3.88 keV with ASCA (Furusho et al. 2001). Unless oth- 
erwise specified, the errors are within the 90% confidence 
region for a single parameter of interest. 

2. Observation and Data Reduction 

One offset and three central regions of the Centaurus 
cluster were observed in December 2005 and 2007 by 
Suzaku. The observed regions range to 23', which cor- 
responds to ~ 0.17 ri 8 o, from the cD galaxy NGC 4696 at 
the center. The CEN45 field corresponds to the subclus- 
ter, the Cen 45 region. The observation logs are given in 
table 1, and an XIS image in 0.5-4.0 keV is shown in fig- 
ure 1. We used only the XIS data in this study. The XIS 
instrument consists of four sets of X-ray CCD (XIS 0, 
1, 2, and 3). XIS 1 is a back-illuminated (BI) sensor, 
while XIS 0, 2, and 3 are front-illuminated (FI) sensors. 
XIS 2, however, has not been available since November 
2006; hence, the CEN 45 observation was carried out only 
by XIS 0, 1, and 3. The XIS was operated in the normal 
clocking mode with standard 5x5 or 3x3 editing modes. 

The analysis was performed with HEAsoft version 6.7 
and XSPEC 12.5.1. We used the standard data selection 




Fig. 1. XIS image (0.5-4 keV) of the Centaurus clus- 
ter. The images obtained with the four pointing ob- 
servations are superposed. The exposure was corrected, 
but the NXB and CXB were not subtracted. The image 
was smoothed by a Gaussian of a =10 pixels. The unit 
of color bar is counts pixel -1 s _1 . The circles indicate 
the annular regions as mentioned in the spectral analysis. 

\ and the CALDB files version 2009-04-02. Because the 
energy resolution slowly degraded after the launch because 
of radiation damage, this effect was included in the redis- 
tribution matrix file generated by the "xisrmfgen" Ftools 
task. We generated an Ancillary Response File (ARF) for 
the spectrum of each annular sky region, which assumed 
r = 25' circle size of the surface brightness profile, which 
is a sum of two /3-models, (3\ =0.57 ,r c =7.'3 and (32 =0.92, 
r Ci2 =0f92 derived from ASCA and ROSAT results (Ikebe 
et al. 1999) by the "xissimarfgen" Ftools task (Ishisaki 
ct al. 2007). We also included the contamination effect 
on the optical blocking filter of the XIS in the ARFs. We 
employed the night Earth database generated by the "xis- 
nxbgen" Ftools task to subtract non- X-ray background 
(NXB). 

3. Spectral Analysis 

3.1. Estimation of the Cosmic X-ray Background and 
Galactic Foreground 

We estimated the cosmic X-ray background (CXB) us- 
ing four Lockman hole observations with Suzaku as shown 
in table 1. Because the observed emissions in our ob- 
servation were dominated by the ICM emission from the 
Centaurus cluster, we could not estimate the CXB level 
in our observation. Thus, we used the Lockman hole data 
for the background estimation. We fitted the spectra 
from the Lockman hole data with an absorbed power- 
law model. The resultant parameters are shown in ta- 
ble 2. The weighted average of photon index T = 1.47, 
and Norm= 9.9 photons cm~ 2 s _1 sr _1 keV -1 at 1 keV 

1 http://www.astro.isas.ac.jp/suzaku/process/v2changes/ crite- 
ria_xis.html 
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Table 1. Suzaku observations of the Centaurus cluster and the background region. 



Target name 


Sequence 


Date 


Coordinates (J2000.0) 


Exposure 




Number 




a (degree) 


S (degree) 


(ks) 


CENTAURUS_CLUSTER 


800014010 


2005-Dec-27 


192 2054 


(1 001 1 1 

—41 3111 


36.3 


CEN CLJJl b SE I 1 


800015010 


2005-Dec-28 


192°2054 


-41°4440 


44.7 


CENCL_0FFSET2 


800016010 


2005-Dec-29 


192°2054 


-41°1780 


43.0 


CEN45 


802008010 


2007-Dec-24 


192°5119 


-41°3865 


58.2 


LOCKMANHOLE 


100046010 


2005-Nov-14 


163°4063 


57°6108 


77.0 


LOCKMANHOLE 


101002010 


2006-May-17 


162°9366 


57°2557 


80.4 


LOCKMANHOLE 


102018010 


2007-May-03 


162°9257 


57°2581 


96.1 


LOCKMANHOLE 


103009010 


2008-May-18 


162°9369 


57°2546 


83.4 



Table 2. Results of spectral fits of the Lockman hole region for the estimation of the CXB component and the Centaurus offset 
regions for the Galactic components. 



CXB 










Scq. Number 


F 


Norm* 


X 2 /d.o.f. 




100046010 
101002010 


1 co+0.03 

1.00_Q Q 3 

i 41+0.04 


14 2+ - 6 

q q+o.4 
y.o_ 4 


1120/945 
831/765 




102018010 
103009010 


1 43+0.05 

1 4Q+0.05 


o O+0.4 
8.8_0.4 

q n+0-5 


662/617 
670/529 




Galactic 










OFF1 & OFF2 










(keV) 


Norm MWR 

(xlO- 3 ) 


fcT L HB 

(keV) 


Norm WB 

(xHT 3 ) 


X 2 /d.o.f 


W +t) m 

u - io -0.03 


O 1+1. V 


0.1 (fixed) 


0+° 
u - u -o.o 


535/561 



* The units of photons cm 2 s 1 sr 1 keV 1 at 1 keV. 

t Normalization of the apec component divided by the solid angle, f2 u , assumed 
in the uniform-sky ARF calculation (20' radius), Norm = J n c nndV / [4tt (1 + 
z) 2 Dj^]/£l { xl0 — 14 cm -5 4007T arcmin - 2 , where D\ is the angular distance to the 
source. 



were used. 

The Galactic emissions mainly arise from the local hot 
bubble (LHB) and the Milky Way halo (MWH). We used 
a two-temperature model {apec; Smith et al. 2001) for 
these emissions. The version 1.3.1 of APEC code is used 
through this paper. For the outermost region (15'-23') 
in the offset 1 and 2 regions, we fitted the spectra with 
the model formula apecLHB + wabs * (apecMWH + power- 
lawcxB + apecicM) as described in Sato et al. (2007a) and 
Komiyama et al. (2009). Here we fixed the LHB temper- 
ature at 0.1 keV that was typical value of XMM result 
(Lumb et al. 2002), while the MWH temperature was a 
free parameter. The normalizations of the LHB and MWH 
were allowed to vary. We also assumed a solar abundance 
and zero redshift for the LHB and MWH components. 
The resultant parameters are shown in table 2. The de- 
rived values were consistent with those of Yoshino et al. 
(2009). 

Using these parameters, we generated mock spectra in- 
cluding the CXB, the LHB, and the MWH components 
as the background. In the spectral analysis of the in- 
ner (0'-15') regions, we subtracted the mock spectra as 
the background from the observed spectra for each sensor. 



The CXB and Galactic components in 15'-23' region were 
included as models in the spectral fits. As a result, the 
temperature and abundances in the ICM did not change 
within the statistical errors when we examined the sys- 
tematic errors of our results by changing the background 
normalization by ±10% as mentioned in subsection 4.5. 

3.2. ICM Component 

We extracted spectra from six annular regions of 0'- 
2', 2'~4', 4'-6', 6'-9', 9'-15', and 15'-23', centered on 
the peak of X-ray intensity at (RA, Dec) = (192°2058, 
-41°3042) (Ota et al. 2007). The annular spectra are 
shown in figure 2. Each spectrum was binned carefully to 
observe details in metal lines, especially below ~ 1 keV. 
Each spectral bin contained 50 or more counts. We fitted 
the spectra in 0.5-9.0 keV at XIS 0, 2, 3 and 0.4-9.0 keV 
at XIS 1. We excluded the narrow energy band around 
the Si-K edge (1.82-1.84 keV) because of an anomalous 
response. 

We assumed a single- or two-temperature (hereafter, IT 
or 2T, respectively) model (vapec; Smith et al. 2001) and 
fitted the spectra of all detectors simultaneously for each 
region. We also examined a three-temperature (3T) model 
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Energy (keV) Energy (keV) Energy (keV ) 



Fig. 2. The panels show the observed spectra for the annular regions of the Centaurus cluster. The data are plotted by black, 
red, and blue crosses for XIS 1 spectra of the CENTAURUS.CLUSTER, CENCL.OFFSET1, and CENCL.OFFSET2 region in 
0'-9', while CEN45, CENCL.OFFSETl, and CENCL.0FFSET2 in 9'-23', respectively. The estimated NXB, CXB, and Galactic 
components were subtracted from the spectra in the 0'-15' region. In 15'-23', the CXB component is a fixed Lockmanhole value 
and Galactic components are included as the models. The yellow lines show the best-fit model for XIS 1. The XIS 1 spectra of the 
ICM components are shown in magenta, orange, and light blue lines. The energy range around the Si-K edge (1.82-1.84 keV) was 
ignored in the spectral fits. The lower panels show the fit residuals in units of u. 



at the innermost region within 2', because the Chandra 
result suggested that at least a 3T model was needed to 
represent the spectra in the central region (Sanders et al. 
2006). 

The metal abundances of He, C, N, and Al were fixed to 
a solar value, while the other nine elements were allowed 
to vary. Note that Ne-K and Ni-L lines could not be re- 
solved from the Fe-L lines. Therefore, we did not report 
Ne and Ni abundances in this paper. The Galactic ab- 
sorption for neutral hydrogen TVh was also allowed to vary 
in the spectral fits. This is because the galactic latitude 
of the Centaurus, 21°7, is relatively low, and there may 
be systematic uncertainties in the adopted Galactic ab- 
sorption which is a weighted average of those of positions 
within 1 degree from the cD galaxy. In addition, a sys- 
tematic uncertainty in the thickness of the contaminant 
on the optical filter of the XIS also causes a systematic 
uncertainty in the Nr. 

4. Results 

4-1. Fitting Results 

Tabic 3 and figures 3, 4, and 5 summarize the resultant 
parameters of the fits. The data-to-model ratios of the fits 
are also shown in figure 3. Within 2', or r < 0.015 7*180, the 
reduced x 2 with the IT model is much larger than unity, 
and the error ranges of parameters are not represented. 
The 2T model significantly improved the resultant \ 2 of 
the fit compared to that with the IT model. The fit with 
the 3T model is even better than that with the 2T model, 




Energy (keV) 

Fig. 3. The data-to-model ratios of the fits of the XIS 1 spec- 
tra with the IT (black), 2T (red), and 3T (blue) models. 

particularly around 0.8 keV (figure 3). At 2'-4', or 0.015- 
0.03 riso, the 2T model also improved the fits around 0.9- 
1.1 keV. Beyond 4', or r > 0.03 n 8 o, with F-tcsts, the 2T 
model is significantly better than the IT model, consid- 
ering only the statistical errors, and the data-model ratio 
improved by several percent around 1.1 keV. However, 
considering possible systematic uncertainties in the Fe- 
L modeling, the IT model fits still represent the spectra 
fairly well. 



Suzaku observations of the Centaurus cluster 
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Fig. 4. (a) Radial temperature profile derived from the spectral fits with the IT (black cross), 2T (red diamond), and 3T (blue 
diamond) models, (b) Ratios of the normalization of the hotter to the cooler ICM component obtained from the spectral fits with 
the 2T (red diamond) and 3T (blue diamond) models. 



4-2. Temperature profile of the ICM 

The radial temperature profile and the ratios of the nor- 
malizations of the hotter and cooler ICM components from 
the cluster center to r ~ 0.17 rigo are shown in figure 4 and 
table 3. The ICM temperatures of the three components 
at the central region, within 0.015 rigo ~ 2', were ~ 3.8, 
~ 1.7, and ~ 0.8 keV, respectively, and the two higher 
temperatures with the 3T model appear to be similar to 
the hot and cool ICM temperatures with the 2T model in 
the 2'-4' (0.015 < r < 0.03 n 80 ) region. Our results for the 
3T model are consistent with the previous Chandra result 
(Sanders et al. 2006). For the outer region (r > 0.015 rigo) 
the ICM temperature of the hotter component with the 
2T model slightly increased to ~ 4.6 keV in the outermost 
region, while that of the cooler component was almost con- 
stant at 1.7 keV. The resultant temperatures with the IT 
model were consistent with the previous XMM-Newton 
result with the IT model (Matsushita et al. 2007b). A 
comparison of the normalization ratios of the hotter and 
cooler components for the 2T and 3T models is shown in 
figure 4. The ratios reached a peak at ~ 0.04 rigo and 
then decreased to approximately one fourth of its peak 
value in the outermost region. 

4.3. Abundance Profiles of ICM 

Radial abundance profiles of O, Mg, Si, S, Ar, Ca, and 
Fe up to ~ 0.17 rigo are summarized in table 3 and fig- 
ure 5. All the abundances lie over a solar abundance at 
the central region, and decrease toward the outer region. 

Within the cool core region (r < 0.045 rigo ~ 6') the 
abundances sharply increased toward the center. The 
Fe abundance reached a maximum value of ~1.83 solar 
within 0.015 n 80 with the 3T model. The Si, S, Ar, and 
Ca abundances showed the same behavior as the Fe pro- 
file. The O and Mg abundances were approximately 1 
solar within 0.015 rigo. 

For the outer region (r > 0.045rig ), the IT and 2T 
model gave similar O, Mg, Si, S, Ca, and Fe abundances, 



although the 2T model gave slightly higher abundances of 
Ar. All the abundances were nearly constant at ~0.5 so- 
lar. Fe decreased to approximately one fourth of the cen- 
tral value in the outer region. The ratio of O and Mg to Fe 
also decreased to approximately half of the central value. 
In the outermost region (0.11<r<0.17 rigo), the fits gave 
only the upper limit of the O abundance with 90% error. 
It is the first time to derive Mg abundance to 0.17ri 8u . 
The derived abundances are consistent with those from 
XMM and Chandra data (Sanders et al. 2006; Matsushita 
et al. 2007b), but the error bars are small particularly 
outside the cool core. 

4-4- Profiles of ICM Abundance Ratios 

The radial profiles of the abundance ratios of O, Mg, 
Si, S, Ar, and Ca to Fe abundance are shown in figure 
6. Excluding the central region, the IT and 2T mod- 
els gave similar abundance ratios of O/Fe, Mg/Fe, Si/Fe, 
and S/Fe, whereas the 2T model yielded higher Ar/Fe and 
Ca/Fe ratios by several tens of percent. The abundance 
ratios of O/Fe and Mg/Fe are ~0.5 solar ratio within 0.03 
rigo and ~0.7-l solar ratio beyond the radius. The other 
abundance ratios of Si/Fc, S/Fc, Ar/Fe, and Ca/Fe are 
~1 solar ratio in the entire region. The weighted averages 
of the abundance ratios are calculated for three radial re- 
gions and are shown in table 4. 

4-5. Uncertainties 

In order to estimate systematic errors, we varied the 
normalizations of the NXB spectra and the mock spec- 
tra including the CXB, LHB, and MWH by ±10% in the 
spectral fits. The systematic errors due to background 
estimation were almost negligible, and the resultant tem- 
peratures and abundances did not change within the sta- 
tistical errors by changing the background level. We also 
note that Ne abundance is not reliably determined be- 
cause of an overlap with the strong and complex Fe-L 
line emissions. However, these abundances were allowed 
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Table 3. Summary of the resultant parameters of fits to each annular spectrum of the Centaurus cluster 



[Vol. , 



r 










A H 


Norm coo \i 


Norm coo \2 


xVd.o.f. 


(arcmin) 




(keV) 


(keV) 


(keV) 


(10 22 cm- 2 ) 


1 Norrrihot 


1 Norm hot 




0'-2' 


3T 


q 7C+U.U4 
03 


1.70+Hi 


0.79±r l 


n i q 7 +o.oo2 

U.IO 1 Q QQ2 


i oc+u.uy 

-L.OU Q Q2 


0.097™ 


8589/7234 


0'-2' 


2T 


2.41±g : gi 


0.99 + °'°i 




r, 1 q 9 + 0.001 


0.071+Hur 

— U.UU1 




9405/7236 


0'-2' 


IT 


2.11 






0.128 






15904/7238 


2'-4' 


2T 


o 70+0.03 
°' ' ° — o 03 


1 69+ 002 

—0.02 




n i 97+O.OO2 
— 001 


0.43±g;°i 




8505/7490 


2'-4' 


IT 


2 90+° :o1 






n -i 99+O.OO2 

U.±Z^_Q QQ2 






9370/7492 


4'-6' 


2T 


4 41 +0 06 
-Q'PJ 


9 -i 9+0.30 
■ a - L,£ -0.19 




n 1 1 s+0.002 


47+0 02 

u -^'-0.02 




6794/6423 


4'-6' 


IT 


^• 48 -0.03 






n 1 1 c+0.002 
u - iio -0.002 






6981/6425 


6'-9' 


2T 


*- u '-0.03 


1 66+ 20 




120+ 002 

U.1ZU_Q QQ2 


o is +0 - 01 

u - lo -0.02 




7475/6877 


6'-9' 


IT 


q R q+0.03 
°- uo -0.03 






n i i 7+O.OO2 
' -0.002 






7587/6879 


9'-15' 


2T 


4 iy+0-14 


i cq+0.13 




0.109+.K 


12+ - 05 

u - iz -0.05 




7470/6818 


9'-15' 


IT 


O 7Q + 0.03 

°-' y -0.03 






o.io6±8:gg 






7562/6820 


15'-23' 


2T 


4 64+ 31 

^• u ^-0.27 


1 61+ - 24 

1 - ui -0.23 




119+0013 
u.iiz_ 012 


n 99+O.O8 
u -^-0.08 




2386/2254 


15'-23' 


IT 


q oq+0.08 






o nqq+ 0011 






2436/2256 



r 







Mg 


Si 


S 


Ar 


Ca 


Fe 


(arcmin) 




(solar) 


(solar) 


(solar) 


(solar) 


(solar) 


(solar) 


(solar) 


0'-2' 
0'~2' 
0'-2' 


3T 
2T 
IT 


1 00 +(U)8 

i - uu -0.08 
1 14+0-07 
i - 1 ^-0.07 

0.38 


1 - uz -0.05 

1 04+ 006 
0.43 


i cq+0.03 
1 - do -0.03 

0.94 


i cc+0.04 
1 - oo -0.04 
1 4Q+0.04 
1 '^ a -0.04 

0.94 


i -'^-0.11 

1 40+ - 10 
i -^ u -o.io 

0.85 


9 i 7+O.I6 
z - i ' -0.16 
i 70+O.I2 
i -' 8 -0.13 

1.58 


1 oq+0.02 

i - oo -0.02 

2 09+ 05 
z - ua -o.oi 
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to vary freely during the spectral fits. 
5. Discussion 

5.1. ICM Abundance Pattern and Contributions from 
SN la and SN II 

Figure 7 summarizes the abundance pattern of the 
O/Fe, Mg/Fe, Si/Fe, S/Fe, Ar/Fe, and Ca/Fe within 0.05 
ri8o ) 0.05-0.1 riso, and beyond 0.1 rigo of the Centaurus, 
AWM 7, Abell 262, and the Perseus clusters with those of 
the four groups of galaxies and the Fornax cluster. Within 
0.1 riso, the Centaurus cluster has the smallest error bars 
for O/Fe and Mg/Fe ratios compared to the other clusters 
and groups of galaxies. 

Outside the cool core, beyond 0.05 r^o, the O/Fc, 
Mg/Fe, Si/Fe, and S/Fe ratios of these groups and clus- 
ters of galaxies are close to the solar ratio. The scatter in 
the Mg/Fe ratio might be caused by uncertainties in the 
Fe-L atomic data, because the Mg-K lines are surrounded 
by the Fe-L lines. The Centaurus cluster and the Perseus 
cluster have similar Ar/Fe and Ca/Fe ratios. There is no 
significant difference in the abundance pattern between 
the clusters and groups of galaxies. 



Within 0.05 risoi or the cool core region, Mg/Fe ratio 
of the Centaurus cluster shows the lowest value of 0.5 in 
units of the solar ratio. Mg/Fe ratio of the other systems 
scatters around unity, i.e., the solar ratio. O/Fe ratio 
of the Centaurus cluster is the same as those of other 
clusters and groups of galaxies, although error bars of the 
other clusters and groups are fairly large. In contrast, the 
Si/Fe and S/Fe ratios of the Centaurus clusters and other 
systems agree well. 

5.2. Mass-to-Light Ratios 

The metal MLR is a useful measure for the study of 
ICM chemical evolution because metals in the ICM are 
all synthesized in galaxies. We calculated B-band and 
K-band luminosity profiles to estimate the MLRs. The 
K-band luminosity of a galaxy correlates more than the 
B-band with stellar mass. The B-band luminosities were 
used to compare previous studies. 

To derive a B-band luminosity profile, we integrated 
the luminosities of the member galaxies from Jerjen and 
Dressier (1997), deprojected the profile assuming spher- 
ical symmetry, and derived a three-dimensional profile 
of Lb- Jerjen and Dressier (1997) optically studied the 
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Fig. 5. (a) Radial neutral hydrogen column density (Nn) profile derived from the spectra. Black dashed lines show Galactic TVjj 
of 8.56 X fO 20 cm -2 from Kalberla et al. (2005). Black crosses, red and blue diamonds correspond to the results with the IT, 2T, 
and 3T models, respectively, (b)-(h) Radial abundance profiles are plotted similar to that in (a). 



Centaurus cluster and identified 296 member galaxies 
whose apparent magnitudes are brighter than 21.5 in the 
B-band. The brightest cluster galaxy NGC 4696 has an 
apparent magnitude tob = 11.3 or logLa/LB,© = H-2 
using 44.9 Mpc and the foreground Galactic extinction 
A B = 0.492 (Schlegel et al. 1998) from NASA/IPAC 
Extragalactic Database (NED). 

To derive a K-band luminosity profile, we used the 
luminosity of galaxy data in a box of 10x10 deg 2 cen- 
tered at the peak of X-ray intensity from the Two Micron 
All Sky Survey (2MASS) 2 . The brightest cluster galaxy 
NGC 4696 has an apparent magnitude uik = 7.298 or log 
Lk/Lk,© = 11.7 using the foreground Galactic extinction 
A K = 0.042 (Schlegel et al. 1998) from NED. Figure 8 
shows the galaxies detected with 2MASS within ~0.2 rigo 
from the cD galaxy. The average surface brightness in 
the region of 140' < r < 300' (ri 80 ~ 140') was subtracted 
as the background. Then, we calculated the deprojected 
three-dimensional profile of Lk from the brightness profile 



http : / / www .ipac.caltech.edu/2mass/ 



of L K . 

Figure 9 shows the integrated radial profiles of Lb and 
Lk- Within 0.1 ri 80 , the cD galaxy dominates the Lb and 
Lk profiles. Uncertainties in the three-dimensional posi- 
tion of bright galaxies except the cD galaxy cause system- 
atic uncertainties in the three-dimensional profile of Lb 
and Lk- At 15' (~0.1 JTgo) offset from the cluster center, 
a luminous elliptical galaxy NGC 4709 exists in the middle 
of the CEN 45 field of view. L K of this galaxy is 40% of 
that of the cD galaxy and causes a sudden increase in the 
integrated Lb and Lk profiles. However, the heliocentric 
radial velocity of NGC 4709 is more than 1 500 km / s higher 
than that of the cD galaxy NGC 4696. Distance measure- 
ments with the surface brightness fluctuation and a glob- 
ular cluster luminosity function indicate that NGC 4709 
is located in front of the Centaurus cluster (Mieske et al. 
2005). Therefore, excluding NGC 4709, we deprojected 
the brightness profiles of Lb and Lk and derived the inte- 
grated profiles, which are plotted as dashed lines in figure 
9. The difference between the solid and dashed lines in 
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Fig. 6. Radial profiles of the abundance ratios of O, Mg, Si, S, and Ar to Fe in solar units (Lodders 2003). The colors are the same 
as in figure 5. 



the figure shows the systematic uncertainty caused by the 
uncertainty in the position of NGC 4709. 

We used the radial profile of gas mass derived from 
ASCA and ROSAT observations up to 22' (Ikebe et al. 
1999) and calculated an integrated mass profile of O, Fe, 
and Mg (figure 9). We then derived integrated profiles of 
MLRs for O, Mg, and Fe (OMLR, MMLR and IMLR) us- 
ing B-band and K-band luminosities (figure 10 and table 
5). The error bars of the MLRs include only the abun- 
dance errors. The IMLR, OMLR, and MMLR profiles 
increase with radius to ~ 0.17 riso, although the pro- 
files have a flatter slope in the 0.1-0.17 rigo region where 
NGC 4709 is probably located in front of the Centaurus 
cluster. Excluding the galaxy, the MLR profiles become a 
smoother slope to ^0.17 r 180 . 

Figure 10 compares the integrated IMLR, OMLR, and 
MMLR profiles of the Centaurus cluster and two other 
clusters Abell 262 and AWM 7 and two groups of galaxies 
NGC 5044 and NGC 1550 with dominant elliptical galax- 
ies in their centers. Within 0.1 r 180 , where the central 
galaxies dominate the integrated Lk profiles, the IMLR, 
OMLR, and MMLR of the Centaurus cluster tends to 
be higher than those of the other systems by a factor 
of 1.5 ~ 2. In contrast, at 0.1-0.17 rigo, MLRs of the 
Centaurus cluster become consistent with those of the 
other systems. 

5.3. Metal Enrichment Histories Outside the Cool Core 

Outside the cool core, from 0.05 ri 80 to ^0.17 ri 80 , 
the observed Fe abundance profile of the Centaurus clus- 
ter agrees well with those of the weighted average of 
the nearby cD clusters observed with XMM (Matsushita 



2011). At this radial range, the IMLR, OMLR, and 
MMLR of the Centaurus cluster are also similar to those of 
AWM 7, Abell 262, and NGC 1550 observed with Suzaku, 
considering the difference in the integrated values within 
0.05 rise) . The abundance patterns of O/Mg/Si/S/Fe of 
these systems are almost similar, although the Mg/Fe ra- 
tio of the AWM 7 cluster at 0.05-0.1 r lso tends to higher. 
These results indicate that outside the cool core regions 
of clusters, the metals in the ICM may be universal in 
clusters of galaxies, and clusters of galaxies may have uni- 
versal metal enrichment histories. 

The observed abundance pattern of 
O/Mg/Si/S/Ar/Ca/Fe beyond 0.05 r 180 is between 
those of SN II and SN la from nucleosynthesis models 
(figure 7). Therefore, both SN la and SN II products have 
been mixed into the ICM. The observed solar abundance 
pattern from O to Fe indicates that a large quantity Fe 
was synthesized by SN la, and the number ratio of SNe II 
to SNe la is estimated to be 3^4 (Sato et al. 2007b; Sato 
et al. 2008). The similarity of abundance pattern of these 
clusters and groups of galaxies outside the cool cores 
indicates that the contributions of the two types of SN to 
the metals in the ICM are similar. 

Renzini (2005) found that the OMLR of the ICM is 
a very sensitive function of the slope of the initial mass 
function (IMF). Adopting a Salpeter IMF, the expected 
value is ~0.1 Mq/Lb,©, and an increase in the slope of 
the IMF decreases the value of the OMLR. In contrast, 
the integrated OMLR at 0.17 ri§o of the Centaurus clus- 
ter, excluding NGC 4709, is 0.01-0.04 M Q /L B ,©. We ex- 
pect the OMLR to become 0.02-0.04 M q /Lb,o, assuming 
that the ratio of O and Mg is constant, since within the 
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Fig. 7. Weighted averages of the abundance ratios of O, Mg, Si, S, Ar and Ca to Fe for the three regions (a) r < 0.05 rigo, (b) 
0.05 < r < 0.1 ri8o, and (c) r > 0.1 riso- We used the abundance ratio derived from the 3T model within 0.015 riso and the 2T model 
for the outer regions. The abundance ratios of other clusters, AWM 7 (Sato et al. 2008), Abcll 262 (Sato et al. 2009b), the Perseus 
cluster (Tamura et al. 2009), NGC 5044 (Komiyama et al. 2009), NGC 1550 (Sato et al. 2010), HCG 62 (Tokoi et al. 2008), NGC 507 
(Sato et al. 2009a), and the Fornax cluster (Matsushita et al. 2007a) are also plotted in the panels. SN II yields by Nomoto et al. 
(2006) refer to an average of the Salpetcr initial mass function of stellar masses from 10 to 50 M® with a progenitor metallicity of 
Z = 0.02. The SN la yields were taken from the W7 model (Iwamoto et al. 1999). 



Table 4. Weighted averages of the abundance ratios in solar units. 
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radius, the abundance ratio of the O and Mg is consis- 
tent with a constant. This value is smaller than the ex- 
pected value obtained with the Salpeter IMF by a factor 
of 3-5. Considering that the observed MLRs of AWM 7 
and Abcll 262 increase with radius, the observed values 
give lower limits of the OMLRs in the entire clusters. 
Therefore, to study the slope of the IMF in clusters of 
galaxies we need measurements of MLRs at outer regions 
of these systems. 



5.4. Metal Enrichment from the cD Galaxy 

The metals in the cool core of clusters are a mixture 
of those in the ICM and the supply from the cD galaxy, 
which contains those synthesized by SN la and those com- 
ing from stars through stellar mass loss. The metal abun- 
dances in the ISM of elliptical galaxies provide impor- 
tant information regarding current metal supply into the 
ICM. In these galaxies, the O and Mg abundances in the 
ISM should be equal to those in mass-losing stars, because 
these elements are not substantially synthesized by SN la. 
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Fig. 8. Identified objects from the Two Micron All Sky 
Survey (2MASS) in K-band (black open circle), the red 
and blue crosses denote the coordinate of the cD galaxy 
(NGC 4696, Ota et al. 2007) and NGC 4709, respectively. 




Fig. 9. Integrated mass of the hot gas (magenta), O (green), 
Mg (red) and Fe (black). Integrated K-band luminosities of 
galaxies Lk (blue) and B-band luminosities of galaxies Lk 
(light blue) are also plotted in the panel. Integrated lumi- 
nosities excluding NGC 4709 are shown with dashed lines. 



Fe abundance of the ISM is a sum of stellar metallicity and 
SN la contribution, which is proportional to -Mg^#gN/a* 
(see Matsushita et al. 2003 for details). Here, is Fe 

mass synthesized in one SN la, £?sn is the SN la rate, and 
a* is the stellar mass loss rate. 

Figure 11 compares the central abundance pattern of 
the Ccntaurus cluster within 2' with those of early-type 
galaxies observed with Suzaku. O/Fe and Mg/Fe ratios 
of the Centaurus cluster, 0.5-0.6 solar ratio, are system- 
atically smaller than those of the galaxies at 0.8-1.0 solar 
ratio. The stellar metallicity of NGC 4696, the cD galaxy 



of the Centaurus cluster, calculated from the gradient of 
optical Mg2 indexes within the effective radius is similar 
to those of NGC 720, NGC 4636, and NGC 5044 plotted 
in figure 11 (Kobayashi & Arimoto 1999). Therefore, the 
difference in the abundance ratios may be caused by a dif- 
ference in contribution from SN la, and the metals in the 
center of the Centaurus cluster are not simple accumu- 
lations of hot ISM in ellipticals. Because the abundance 
pattern outside the cool core regions of clusters is similar 
to those of ISM in elliptical galaxies, the metals in the 
center of the Centaurus cluster may not be a simple mix- 
ture of those in the ICM and the current ISM in elliptical 
galaxies. 

The difference between the Centaurus cluster and el- 
liptical galaxies is the accumulation time scale of metals 
from stars. The IMLR using B-band luminosity of ellipti- 
cal galaxies observed with ASCA is 10~ 5 - 1O~ 4 M /L B . Q 
(Makishima et al. 2001), which is a factor of 3-30 smaller 
than that within 2' (0.015 rigo) of the Centaurus cluster. 
The enrichment timescale of excess Fe in the center of the 
Centaurus cluster is at least several Gyr (Bohringcr et al. 
2004), whereas the enrichment time scales of hot ISM in 
elliptical galaxies are smaller than 1 Gyr (Matsushita et 
al. 2000; Matsushita 2001). Therefore, a simple interpre- 
tation is that a longer time scale indicates that the ratio 
of SN la rate to stellar mass loss rate was higher in the 
past. 

Another interpretation is that SN la products in ISM 
in elliptical galaxies are lost to intergalactic space by their 
buoyancy, as discussed in Matsushita et al. (2000), based 
on ASCA observations of early-type galaxies. Tang and 
Wang (2010) simulated the evolution of hot SN la ejecta 
and found that they quickly reach a substantially higher 
outward velocity than the ambient medium. In the center 
of the Centaurus cluster, these ejecta can mix with the 
surrounding medium because of the long enrichment time. 

The central value of Fe abundance of the Centaurus 
cluster is highest among nearby cool core clusters (Sanders 
et al. 2006; De Grandi and Molendi 2009; Matsushita 
2011) within 0.03 rigo, Fe abundance of the Centaurus 
cluster is a factor of two higher than the average of the 
other clusters with cool cores. Within the cool core, the 
slope of Fe abundance of the Centaurus cluster is the 
steepest, while those of the other nearby clusters with the 
cool cores are similar (Matsushita 2011). 

With XMM observations of the Centaurus, Virgo, 
Perseus, and Abell 1795 clusters, Bohringer et al. (2004) 
found that enrichment times for the central 20 kpc of the 
Centaurus cluster is a factor of two greater than those of 
the Virgo and Perseus clusters and a factor of 5 greater 
than that of the Abell 1795 cluster. We found that the 
integrated IMLR of the Centaurus cluster is a factor of 
two higher than that of Abell 262 and AWM 7 clusters. 
Furthermore, the central O/Fe and Mg/Fe ratios of the 
Ccntaurus cluster tend to be smaller than those of the 
other groups and clusters of galaxies (figure 7). These 
results indicate that the central regions of the Centaurus 
cluster contain more SN la yields than the other clusters. 

With Chandra observations, Fabian et al. (2005) found 



Suzaku observations of the Centaurus cluster 

Table 5. Summary of B-band and K-band metal mass-to-light ratios of the Centaurus cluster 
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Fig. 10. (a) Integrated Fe mass-to-light ratio (IMLR), (b) Oxygen mass-to-light ratio (OMLR), (c) Magnesium mass-to-light ratio 
(MMLR) in units of Mq/Lq in the B-band (cross) and K-band (open circle). Dashed line represents the result using the calculated 
luminosity, excluding NGC 4709. For comparison, we plotted the radial profile of cluster of galaxies AWM 7 (Sato et al. 2008, Sato 
et al. in preparation) and Abell 262 (Sato et al. 2009b) and galaxy groups NGC 5044 (Komiyama et al. 2009) and NGC 1550 (Sato 
et al. 2010). 



that the core of the Centaurus cluster is complicated with 
bubbles and filaments, and sloshing motions of the gas 
within the central potential well or the central active 
galactic nuclei may cause the central structure of the ICM. 
Dynamics of the ICM in the cool core, particularly turbu- 
lent motions, have been used to explain the lack of cool 
X-ray emitting gas in the cores of galaxy clusters. The 
turbulent motion in the core of the Centaurus cluster is 
loosely constrained as <1400 km s" 1 from the Suzaku/XIS 
analysis of Fe-line width for the Doppler broadening (Ota 
et al. 2007) and <1100km s~ 4 from the XMM/RGS ob- 
servations (Sanders et al. 2011). Graham et al. (2006) 



derived the limit of 400 km s~ 4 in the central 25 kpc of 
the cluster based on an argument that strong turbulence 
would smear out, via diffusion, the observed abundance 
gradient. A tighter constraint will be attained with the 
future high-resolution spectrometer onboard the ASTRO- 
H (Mitsuda et al. 2010; Takahashi et al. 2010). The larger 
contribution from SN la limits the mixing in the cool core 
in the center of the Centaurus cluster, compared to that 
in the other cool cores. 
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Fig. 11. A Comparison of the abundance ratios in the 
0' < r < 2' region of the Centaurus cluster (black) with the ra- 
tio of SO galaxies (blue) and elliptical galaxies (red). SO galax- 
ies are two-temperature results of NGC 1316 (Konami et al. 
2010) and NGC 4382 (Nagino et al. 2010). Elliptical galaxies 
are NGC 720 (Tawara et al. 2008), NGC 1399 (Matsushita et 
al. 2007a), NGC 1404 (Matsushita et al. 2007a), NGC 4636 
(Hayashi et al. 2009), and NGC 5044 (Komiyama et al. 2009). 

6. Summary and Conclusion 

The abundance pattern of O, Mg, Si, S, Ar, Ca, and 
Fe are derived up to the outermost radius 0.17 ?*iso 
with Suzaku observations. Outside the cool core (r > 
0.07 rigo) all elements have almost the same value of 
0.5 solar and the abundance pattern is approximately 
1 solar ratio. The integrated IMLR is similar to the 
other cluster of galaxies. The similarity of abundances 
and MLRs indicates that clusters of galaxies may have 
universal metal enrichment histories outside the cool 
cores. Within the cool core (r <0.05 rigo) the abundances 
suddenly increase toward the center. Si, S, Ar, Ca, and 
Fe abundances have a value of 1.5 ~ 2 solar, while O and 
Mg abundances have a value of ~1 solar. The highest 
central Fe abundance and IMLR and smaller O/Fe and 
Mg/Fe ratios compared to the other clusters limit the 
mixing in the cool core in the Centaurus cluster. 
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